Quantum discord for general two-qubit states: Analytical progress 
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We present a reliable algorithm to evaluate quantum discord for general two-qubit states, amending and 
extending an approach recently put forward for the subclass of X~states. A closed expression for the discord 
of arbitrary states of two qubits cannot be obtained, as the optimization problem for the conditional entropy 
requires the solution to a pair of transcendental equations in the state parameters. We apply our algorithm to 
run a numerical comparison between quantum discord and an alternative, computable measure of non-classical 
correlations, namely the geometric discord. We identify the extremally non-classically correlated two-qubit 
states according to the (normalized) geometric discord, at fixed value of the conventional quantum discord. The 
latter cannot exceed the square root of the former for systems of two qubits. 

PACS numbers: 03.67.-a, 03.65.Ta 



I. INTRODUCTION 

It is as well theoretically proven as experimentally tested 
that the exploitation of entangled states improves our ability 
to manage information in several ways [T]. The development 
of quantum information theory is definitely due to this funda- 
mental result 1 2]. However, it has been recently discovered 
that even multipartite separable states can play a relevant role 
in performing better-than-classical communication and infor- 
mation protocols l3[-i5|]. In general, then, the usefulness of 
a particular state of a quantum system for such tasks can be 
traced back to the presence of internal correlations among 
parts of the system, i.e., to the knowledge that an observer 
Alice, probing a subsystem A, gains about the state of another 
subsystem B, controlled by another observer Bob, and vice 
versa. 

The quantitative and qualitative evaluation of such correla- 
tions and the proper discrimination of their nature — clas- 
sical versus quantum — stand as open problems. Several 
quantifiers of non-classicality of correlations have been in- 
troduced in Uterature ll6l-IT5l (not to be confused with the 
non-classicality of quantum optical states fT^), but there are 
still neither clear criteria of faithfulness or bona-fide-ness for 
them, nor a well established hierarchy of reliability. In spite 
of that, they are heavily used in current research [IITI^SOII. The 
most popular one is by far the quantum discord llaHl, a mea- 
sure of non-classical correlations which goes beyond entan- 
glement and whose definition has an immediate interpreta- 
tion in information theory: discord equals the difference of 
the total correlations between two subsystems A and B be- 
fore and after a local measurement process is performed on 
one of them. The quantum discord admits at least three other 
operational interpretations, in contexts ranging from thermo- 
dynamics to communication protocols, such as the quantum 
state merging |21]. The evaluation of the quantum discord 
is a hard task from a computational point of view, implying 
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an optimization of the conditional entropy between Alice and 
Bob, S{A\B), over all local (generalized) measurements on 
one party, which is often obtainable by numerical methods 
only. A closed analytical solution is known in the case of ar- 
bitrary two-mode Gaussian states |20|, under the restriction 
of Gaussian local measurements. Narrowing our overview to 
two-qubit states, an analytical expression of discord has been 
derived in particular for the subclass of so-called X-states 
ll24l - l27ll . and a successful attempt to generalize this proce- 
dure has not been advanced yet (to the best of our knowledge), 
apart from an upper bound quantity for the discord defined in 

m- 

The difficulty in calculating quantum discord motivated the 
introduction of alternative measures of non-classical correla- 
tions. In particular, the geometric discord [13] is one such a 
measure, which quantifies the amount of non-classical corre- 
lations of a state in terms of its minimal distance from the set 
of genuinely classical states. The geometric discord involves a 
simpler optimization and is easily computable analytically for 
general two-qubit states. However, its relationship with the 
original quantum discord is not entirely clear at the present 
stage. 

In this work, we present an algorithm to calculate quan- 
tum discord for general two-qubit states. First, we obtain 
an explicit and simplified expression for the conditional en- 
tropy, exploiting the Bloch representation of the density ma- 
trix; then, we employ new variables that allow us to set the 
optimization conditions in a closed form. Finally, we asso- 
ciate them to constraints over the eigenvalues of the statistical 
ensemble obtained after the measurement process. Our ap- 
proach qualifies as the most efficient and reliable way to eval- 
uate quantum discord for arbitrary states of two qubits fyi\ . 

Exploiting our algorithm, we perform a detailed numerical 
exploration of the Hilbert space of two-qubit states to com- 
pare quantum discord and the geometric discord as quantifiers 
of non-classical correlations. We shed light on the relation- 
ship between these two quantities by identifying the states that 
extremize geometric discord at fixed quantum discord (and 
vice versa). We are motivated by the aim of establishing a 
reliable hierarchy of non-classical states based on physically 
and mathematically consistent criteria. We find that, interest- 
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ingly, the quantum discord of a two-qubit state can never ex- 
ceed its (normalized) geometric discord. In analogy with the 
study of maximal entanglement ll32l] . we find that the notion 
of maximal non-classical correlations is measure-dependent: 
Therefore, the feasibility in a specific experimental realization 
will determine which, among the various classes of maximally 
quantumly correlated states, are the most suitable ones to be 
employed for applications. 

The paper is organized as follows. Section provides an 
introduction to quantum discord and geometric discord (here- 
after denoted by D and Dc, respectively), presenting defini- 
tions, main properties and a summary of the results obtained 
in the literature. In Section |III] we pursue an analytical ap- 
proach to the calculation of quantum discord for general two- 
qubit states, eventually recasting the optimization problem 
into a system of two transcendental equations, whose solu- 
tion specifies the local measurement that minimizes the condi- 
tional entropy. Section |IV] concerns the comparison between 
quantum discord and geometric discord by using the results 
of the previous sections; we identify the classes of states with 
maximal and minimal geometric discord at fixed quantum dis- 
cord. Finally, Section |V]recalls the main results of our work 
and suggests further issues worthy of investigation. 



II. BASIC DEFINITIONS 
A. Quantum discord 

One of the lessons we can apprehend from Quantum Me- 
chanics is that the measurement process disturbs the state in 
which a physical system is set. This differs from what happens 
in the classical scenario, hence it is possible to conclude that 
the disturbance induced by a measurement on a state is a good 
evidence of its "quantumness". Now, let us suppose to have 
a bipartite system in a certain state and make a measurement 
on one of its subsystems. We can analyze the nature of the in- 
ternal correlations of the system in such a state by studying in 
which way they are affected by the measurement. Information 
Theory provides the tools to accomplish this task. 

As a starting point, let us consider the case of information 
stored in two classical probability distributions. The quantity 
expressing the total amount of correlations between two ran- 
dom variables X and Y, assuming values {x,), {yj] with proba- 
bility {/?,), [qj], is the mutual information I, defined as 



I{X : Y) = 'H(X) + 'H(Y) - 'H(X, Y), 



(1) 



where "HiX) = - 2, p, log2 pi is the Shannon entropy as- 
sociated to the random variable X; consequently 'K(F) = 
- Yjj Ij log2 qj is the entropy for Y while '7/(X, Y) is the joint 
entropy of X and Y. Following Bayesian rules, we can retrieve 
equivalent formulations, I(X : Y) = I(Y : X) = J(X : Y) = 
J(Y ; X), where 



J{X : Y) = 'H{X)-'H{X\Y), 
J{Y:X) = 'H{Y)-'H{Y\X). 



(2) 
(3) 



Here the conditional entropy is straightforwardly defined as 
'H{X\Y) - "HiX, Y) - l-liY), and represents how much uncer- 
tainty (ignorance) we have on X given the value of Y (and vice 
versa for 'H{Y\X)). 

In the quantum scenario, we consider a bipartite system AB 
described by a density matrix p = pab, and subsystems A, B 
with marginal density matrices Pa,Pb- The mutual informa- 
tion can be used once more to quantify the total correlations 
between A and B. The quantum analogue of the expression 
([T]i straightforwardly reads 



J(A : B) = S(A) + S(B) - S{A, B), 



(4) 



where S{A) - -Tt[pa log2PA] is the von Neumann entropy 
of subsystem A, and equivalently S{B) = -Tr[pB log2PB], 
S{A, B) = -Tr[plog2p]. On the other hand, if we try to define 
a quantum version of J' we have 



J{A : B) = S(A) - S(A\B), 



(5) 



which is an "ambiguous" quantity |0,|7tl, since the conditional 
quantum entropy S{A\B) is clearly depending on which ob- 
servable we have measured on B. Using the same notation of 
ll25ll . we recall that a von Neumann measurement (from now 
on just measurement) on B projects the system into a statisti- 
cal ensemble {pk,pii}, such that 



•Pk = 



(Ia < S> PBk)p(lA <S> Psk) 
Pk 



where 



(6) 



(7) 



Pk = Tr[p(lA®PBk)] 

PBk = VTlkV^' 

Rk = \k}{k\, k^O,l 

V e SU(2). 



We remind that in this case study the use of generalized 
positive-operator-valued-measurements (POVMs) is not re- 
quired, since it has been proven in [[23] that for two-qubit 
states the optimal measurement for the conditional entropy is 
always a projective one. We can say that the amount of truly 
classical correlations is expressed by the mutual information 
obtained adopting the least disturbing measurement ||7|: 

C(A : B) = maX|p„(J-(A : B) (8) 
= SiA) - minip^^j 2 PkS(A\B^p,,^). 

k 

Consequently, the amount of genuinely quantum correlations, 
called quantum discord, is given by la] 

2)(A : B) = I(A : B) - C(A : B) (9) 
= S{B) - S{A, B) + minip,,, ^ pkS(A\BiP,,0- 



We note that quantum discord is not symmetric: 
S(A)-S{A\B) + S{B)-S(B\A); 



(10) 
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performing the measurement on Alice's subsystem rather than 
on Bob's one is perfectly legitimate, but it returns in general 
a different value of discord. See e.g. lfl4l I22I1 for a discus- 
sion about the implications of it. It is immediate to verify that 
not only entangled states, but almost all separable states have 
a non-vanishing quantum discord lfl9ll . i.e. are affected by 
the measurement process, thus exhibiting some pretty quan- 
tum properties. In the case of pure bipartite states, the dis- 
cord reduces to the marginal entropy of one of the two subsys- 
tems and therefore to the canonical measure of entanglement. 
Quantum discord for two-qubit states is normalized to one. 

B. Geometric discord 



Bloch vectors associated to the subsystems A, B, and f,y denote 
the elements of the correlation matrix T. Then, it is shown in 
ifisll that the geometric discord is given by 

Dc(p)^\(\\yy^\\2 + \\T\\l-k), (14) 

with k being the largest eigenvalue of the matrix yy^ + T^T 
(in case of measurement on Alice, one needs to replace y with 
X and T^T with TT^). An alternative formulation for the geo- 
metric discord has been provided in [30]. It is easy to see that 
Dc is not normalized to one: its maximum value is 1/2 for 
two-qubit states, so it is natural to consider 2Dc as a proper 
measure for a comparison with the quantum discord D. 



Recently, it has been argued that the experienced difficulty 
of calculating quantum discord can be coped, for a general 
two-qubit state, with the introduction of its geometrized ver- 
sion, hereby just called geometric discord 1 13]. 
Let us suppose, to be coherent with section UTAl to have a bi- 
partite system AB and to make a measurement on B. As we 
have remarked, almost all (entangled or separable) states are 
disturbed by the measurement; however, there is a subclass 
of states which is invariant and presents zero discord. It is 
the class of the so-called classical-quantum states [4^, whose 
elements have a density matrix of this form 



P 



^ PiPAi ® \i){i\: 



(11) 



where p, is a probability distribution, pA, is the marginal den- 
sity matrix of A and {|;)} is an orthonormal vector set. A 
classical-quantum state is not affected by a measurement on 
B in any case. 

Letting Q be the set of classical-quantum two-qubit states, 
and be a generic element of this set, the geometric discord 
Dc is defined as the distance between the state p and the clos- 
est classical-quantum state. In the original definition ] 13], the 
(squared) Hilbert-Schmidt distance is adopted. Recalling that 
11^112 - Tr[AA^] is the square of the Hilbert-Schmidt norm of 
a matrix A, the geometric discord has been introduced as 



Dc(p) = min^enWp - xWl- 



(12) 



It is possible to obtain an explicit closed expression of Dc 
for two-qubit states. First, one needs to express the 4x4 den- 
sity matrix of a two-qubit state in the so-called Bloch basis 



p = jYjRijO-i^o-j 



(13) 



ij=0 



i=l 

3 3 
+ ^yjhx2 ® cry + ^ tjjO-i (g) CTyj, 



where Rij = Tr[p(cr; igi (Tj)], ctq = 12x2, cr,- (i - 1,2, 3) are 
the Pauli matrices, x - {x, j,y = {y,) are the three-dimensional 



III. QUANTUM DISCORD FOR TWO-QUBIT STATES 
A. General setting 

Even though quantum discord has an apparently simple def- 
inition ]6], the practice reveals that its explicit evaluation is 
hard to accomplish. In this paper we restrict our attention to 
two-qubit states. An analytical algorithm has been proposed 
for the subclass of states with maximally mixed marginals (de- 
scribed by five real parameters) in [24]. Also, an extension to 
states spanned by seven real parameters, called X-states be- 
cause of the peculiar form of their density matrix (with vanish- 
ing elements outside the leading diagonal and the antidiago- 
nal), has been introduced in ll25ll . and amended by [|27ll . Here, 
we attempt to generalize the procedure to the entire class of 
two-qubit states. 

First, we consider that performing local unitary transfor- 
mations we can recast the density matrix for an arbitrary two- 
qubit state, Eq. (fT3l l. in the Bloch normal form 124, . 3 3.1 

p = j(hx4 + ^ fl/CT; » /2x2 + ^ bil2x2 ® cr,- 



^ C/CT; (S> cr,), 



(15) 



that is a density matrix completely defined by nine real param- 
eters arranged in three 3-dimensional column vectors a = {a,), 
b - {bj} and c = {c,). This follows from the fact that local uni- 
tary operations p' - {Ua® Ub)p{Ua ® Ub)^ correspond to left 
and right multiplication of the Bloch matrix R with orthogonal 
matrices ll33ll . 



R' = 



1 
01 



R 



1 
Ob 



(16) 



with Oa,b g S (9(3). It is then straightforward to obtain the 
normal form of Eq. (flSt : one needs to calculate the singu- 
lar value decomposition of the lower diagonal 3 ® 3 block 
T of R, T = OaCO^^, divide Oa and Ob by their respective 
determinant (to make sure they both have determinant -i-l), 
and then apply Eq. ( fTSI ). The density matrix is correspond- 
ingly transformed into the normal form of Eq. (ITSt , where the 
c, are identified with the elements of the diagonal matrix C, 
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i.e., the singular values of T. Every two-qubit state can be 
then transformed in its simplified normal form by means of 
local unitaries (which preserve entanglement and correlations 
in general, by definition), so we can restrict our analysis to 
density matrices of this type without incurring in any loss of 
generaUty. 

Now, we move to calculate the quantum discord £), Eq. (|9]l, 
for generic states in normal form. The marginal entropy S(B) 
and the global entropy S(A, E) are trivial to obtain. The main 
issue regards the optimization involved in the conditional en- 
tropy. 



pression by exploiting the normal form of the density ma- 
trix, Eq. (fTST l. However, we have to remark that in this pic- 
ture there is still an amount of redundancy [[26t]. A projec- 
tive measurement on a two-qubit state can be characterized 
by two independent variables only, identifiable as the angles 
9 and 0, which parametrize a generic single-qubit pure state 
as If/r) - cos 9\Q) + e"^ sin 9\\), and the Bloch sphere of coordi- 
nates [x, y, z] in this way. 



X — 2j - 2cos6sm6cos(f> 
y - 2/i = 2 cos 6* sin sin 
z = 2k -I = 2cos2 6i- 1. 



(21) 



B. Conditional Entropy 

Firstly, we have to write the conditional entropy in an ex- 
plicit form, adopting for simplicity the notation in ll25ll . We 
remarked that a measurement sends the state p into an en- 
semble {pk,Pk] as expressed in Eq. (|6]l. The entropy of the 
ensemble can be written as 

S = 2/?,>S(A|B,f,,)) = ^/,,^(pi) (17) 

k k 

= poSo + PiSi, 

where So,S i are the entropies associated to po,pi- 
The measurement is defined by the quantity Psk and is conse- 
quently parametrized by the elements of the unitary matrix V, 
which we can write in the basis of the Pauli matrices as 



(18) 



We notice that the real vector {vo,v^ = {v,) has norm one. 
Therefore, it is possible to rearrange the four parameters in 
three variables only, for example in this way. 



h — vqVi + V2V3 



= V1V3 - V()V2 

2 2 
= V0-HV3. 



(19) 



Setting the vectors X = {2j,2h,2k - 1) and m-t - {m,-t) - 
{flj ± CiXi}, we have that, after a straightforward calculation, 
the conditional entropy takes the following expression: 

4 1 [\ l-b-Xl \ l-b-Xl 



1 + 



ml 



l-b-X) 

1 - 



log2 1 + 



l-b-X 



l0g2 1 



\m+\ 



m+ 



l+b-X 



log2 1 + 



l+b-X 



(20) 



This result is consistent with the formula provided in the 
Appendix of [.26.1 . but we have reached here a simpler ex- 



It is immediate to verify that the following constraint holds 

k^+h^ + / = k. (22) 

The algorithm originally designed for X-states in ll25ll is 
flawed just in not considering the mutual dependence of h, j, k, 
resulting reliable only for a more restricted class of states 
identified in [27]. 

The above mapping enables us to re-parameterize the con- 
ditional entropy, Eq. ( |20] |. as a function of the azimuthal and 
polar angles 9, (p; we can then write S{h, j, k) - S(6, (p) and 
perform the optimization of S over these two independent 
variables. 

C. Optimization 

Inspired by | |25|[ . we look for symmetries in the expression 
of the conditional entropy. We notice immediately the invari- 
ance under the transformation — > 9 + ji, which, however, can 
be englobed by the following one 




which corresponds to 



9 -> 9 + K/2. 



(23) 



(24) 



We can appreciate this with an example. Let us pick a ran- 
dom state, such as 



0.437 0.126 + 0.197/ 

0.126-0.1971 0.154 

0.0271 +0.0258; -0.0115+0.0187/ 

-0.274 - 0.0997/ -0.03 1 5 - 0. 1 70/ 



0.0271 - 0.0258/ -0.274 + 0.0997/ \ 

-0.0115 -0.0187/ -0.0315 + 0.170/ 

0.0370 0.00219-0.0367/ 
0.00219 + 0.0367/ 0.372 

(25) 

we can operate with local unitaries on it, obtaining a new state 
p (albeit with the same entropies and discord) described by 
the simplified normal form presented in Eq. (fTSt : we then per- 
form a projective measurement on subsystem B, obtaining an 
ensemble whose conditional entropy is plotted in Fig.[T] 

One sees that there are no further apparent symmetries for 
the conditional entropy. Therefore the analysis so far, while 
not being conclusive, allows us just to refine the problem by 
safely letting the optimization of the conditional entropy to be 
restricted to the interval € [0, 7r/2). To determine the min- 
imum of S, we need to calculate its derivatives with respect 
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FIG. 1; (Color online) Example of conditional entropy <S for a ran- 
dom two-qubit state [Eq. (|25H. The angles 6 and (p parametrize 
the measurement: we can appreciate the symmetry properties of .S 
with respect to such variables, expressed by the invariance S{6, <p) = 
S(6 ± n/2, (p). All the quantities plotted are dimensionless. 



to and 0. The dependence on these variables involves loga- 
rithms of non-linear quantities, so we can not expect to solve 
analytically the problem in any case, whatever ingenious vari- 
ables we might choose. However, we can seek to write the 
two constraints in a compact and elegant form. Let us impose 




(26) 



After a bit of algebra, we obtain 

-S = -/:'-r_)log2[l-/:'-r_] 

+ (1 -p-Hr_)log2[l-p-Hr_] 

+ (1 + pH-r+)log2[l + p + r+] 

+ (1 -Hp-r+)log2[l + p-r+] 

_ 4 + 2(-(l-/,)log2[l-p] 

- (I + p)log2[l + p])]. 



(27) 



Now, we set the partial derivatives to zero, 

■ dS _ dS Sp . dS ar+ OS Or- _ r. . 

ae dp de flr+ oe oe ^ ' 

dS _ dS dp ^ dS dr^ ^ aS di _ _ q 
dip dp dtp dr+ dip dr^ dcp 

Defining the following quantities 



det 



yS = det 



det 



/ 


dp 


dp 






de 


d<p 




dn 


dr^ ' 


\ 


de 


dip 


) 


/ 


dp 


dp ^ 






de 


r30 






dr- 


dr- 




\ 


de 


dip ) 




/ 


dr^ 


dr+ \ 






le 


dip 






dr- 


dr- 






de 


dip ) 





(28) 



after some manipulations, we can write the stationarity condi- 



tions in the following form 

■ ilog ril£^l + ilo? r ('+P)('-^-'--) i P - 0- 



llnd r ^-P-''- 1 _ iloij r u-i-wi-/^-'-; -! a 
y 4 '"S2t J 2^'^&2\-{\-p){\+p-r^)ia+p- 



(29) 



0. 



We see immediately that this system can be further simpli- 
fied to 



= 0; 



(30) 



ilod r '~^~''- l - ilosj r u+wi-p-'--; -! a 
\ 4 ^"e2t J 2 ^"62t(i_p)(i+p_,.^)JQ.+/j. 



+P+Y 



0. 



We can still express these equations as relations among the 
eigenvalues of the ensemble {pk,Pk]- Calling /l^, /ly the eigen- 
values of po and the eigenvalues of p\, we have 



2\ l-p 



1 



(31) 



After some straightforward algebra, one can show that the 
vanishing of the derivatives of S occurs when the following 
constraints are satisfied 



(32) 



A- 



These two transcendental equations can be solved numer- 
ically. They represent the most compact formulation to date 
for the problem of calculating the quantum discord of arbi- 
trary two-qubit states. Let us call i, the solutions obtained, 
coiTesponding to values {6*,, 0,). In order to establish if they 
represent minima of S, we adopt the conventional method and 
evaluate the signature of the Hessian matrix H at the points 
{6i, (pi), and, in case of det// - 0, we study the sign of the 
functions 5; - S(6, (p)-S(6i, 0,). Naming {0mj, <pmj] the angles 
such that H is positive definite or (5, > 0, we clearly have that 
the absolute minimum of the conditional entropy is defined as 



The quantum discord for generic two-qubit states of the form 
Eq. ( fTSl ) finally reads 

D{A : B) = S{B) - >S(A, B) + min|e,„^,^,„^|>S(0„„, 4>„,j). (33) 



IV. COMPARISON BETWEEN QUANTUM DISCORD AND 
GEOMETRIC DISCORD 



We can use our results to compare the quantum discord D 
[6] with the geometric discord Dc \ 1 3i1 for general two-qubit 
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0.2 0.4 0.6 0.8 



FIG. 2: (Color online) Comparison between normalized geometric 
discord IDg and quantum discord D for 10* randomly generated 
general two-qubit states. The dashed line is obtained by taking the 
equality sign in Ineq. ( |35t . Refer to the main text for details of the 
other boundary curves. All the quantities plotted are dimensionless. 



states fi^. We have generated up to 10'' random general two- 
qubit states. After transforming each of them into the normal 
form of Eq. ( fTSl l. we have calculated their quantum discord D, 
as numerically obtained from the algorithm of Section|III] and 
their normalized geometric discord 2Dg- The latter admits the 
following explicit analytic expression for states p in normal 
form, derived from Eq. (fT4t . 



Dc{p)^-{\\bb'\\2 + 



\cc II2 ■ 



(34) 



where k is the largest eigenvalue of the matrix blr + c(F . 

The results are shown in Fig. |2] We notice that physical 
states of two qubits fill up a two-dimensional area in the space 
of the two non-classicality measures, meaning that the two 
impose inequivalent orderings on the set of mixed two-qubit 
quantum states; this is reminiscent of the case of entanglement 
measures (see e.g. fy2\), and a similar feature has been re- 
ported concerning the comparison between discord and other 
non-classicality indicators |14]. Nevertheless, at fixed quan- 
tum discord, the geometric discord admits exact lower and 
upper bounds (and vice versa). We have identified them nu- 
merically. 

Firsts of all, we have observed that the following hierarchi- 
cal relationship holds for arbitrary two-qubit states. 



2Dg > 2) 



(35) 



In other words, the quantum discord for any two-qubit state 
can never exceed the square root of its (normalized) distance 
from the set of classical-quantum states. The corresponding 
boundary curve is plotted as a dashed (gray) line in Fig. |2] 
However, we see that such a bound is tight only in the region 
< 2) < 1/3, in which it coincides with what we refer to 
as branch (i) [see below], while it is not attainable for higher 
degrees of non-classical correlations. The actual, tight lower 
bound in the whole {£), 2Dc] plane accommodates states with 



minimal geometric discord at fixed quantum discord, or equiv- 
alently maximal quantum discord at fixed geometric discord: 
Such extremal states are constituted by the union of four dif- 
ferent families, which sit on branches (i)-(iv) in Fig.|2] 

(i) (Green online) This branch is filled by so-called a states 



Pa - 



' - 

2 ^ 

— 

^ 2 














— 

2 ^ 



0<q;< 1/3, 



(36) 



for which 'Dip a) 
Ineq. ( l35T l. 



a and IDcipa) - cr, thus saturating 



(ii) (Blue online) This small branch is filled by a subclass 
of the two-parameter family 



({\-a)l2 r/2 

fl 



r/2 (1 -fl)/2 j 

5 



(37) 



- < a< 
3 " "14 



V4a - 3fl2 - 1 , 



with r e |^ V4a - - 1, given by the solution 

2rtanh"'( Va-+r-) 



to 



+ ln(-a - r + \) - ln(-fl + r -1- 1) + 

2tanh '(r) = 0. The geometric discord of these states 
is simply 2Dc(p,) - a^, while their quantum discord 
is calculated in |28]. We highlight the presence of the 
"pimple" at the joint between branches (i) and (ii), a re- 
curring feature in the profile of extremal states involv- 
ing quantum discord 1 14, l28ll29ll . 



(iii) (Red online) This branch accommodates asymmetric 
X-states of the form 



Vfl - fl^ - ac 

c 



acOO I — a — c 

2^-- 



(38) 



1 -2c-H 



2c 



0<g< 1, 



with 



c € 



1- 



1 . 



otherwise; 



solution to 8(1 - 2c)c^ tanh"' ( y8(c -\)c-2g + 'i) - 
4c2^8(c- l)c-2g-H3tanh"'(l - 2c) + 
2V8(c-c2)-2g + 3(2c2 +g-l) tanh-' (l£^^ffl£_i) 



0. 



For these states. 



J and 

Dipg) = - ln(-4c(fl + c - 1)) - 

2 V4c(flH-c- l)-i- 1 tanh"' ( ^fAc{aT7^^\)TVj 
21n(l - fl) -H 4fltanh"'(l - 2fl) + 21n(2 - 2c) - 
4ctanh-'(l -2c) . 
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(iv) (Black online) The top-right-most branch accommo- 
dates just pure states pp - for which the 
discord equals the marginal von Neumann entropy, 
Dipp) =S(pa) = -/7log2/:'-(l-/7)log2(l-/7), andthe 
geometric discord equals the marginal linear entropy, 
2Dg(Pp) = 2(1 - Ti-pl) = 4detpA = 4p(l - p), where 
we have denoted the eigenvalues of the reduced density 
matrix pA by {p, 1 - p). 

On the other hand, the upper boundary (v) in Fig. |2] despite 
being single-branched, is more involved and we are unable to 
provide a tractable parametrization of the states that saturate 
it. They can be sought among symmetric X-states of full rank, 
but with the two biggest eigenvalues dominating the other two. 
The extremal curve has been obtained as the result of exten- 
sive numerical optimization, in which the parameter space has 
been finely sliced in discrete intervals of nearly constant dis- 
cord, and for each interval the datapoint corresponding to the 
random state with the maximum geometric discord has been 
selected. Joining all such extremal states we have obtained the 
smooth (Magenta online) line of Fig.|2] 

The two measures D and 2Dg correctly coincide on 
classical-quantum states Eq. ( fTTT i. where both vanish, and on 
maximally entangled Bell states, where both reach unity. 

V. CONCLUSIONS 

We have presented a reliable and efifective algorithm for the 
evaluation of the quantum discord D of general two-qubit 



states. We have simplified the optimization involved in cal- 
culating the conditional entropy, by removing the redundant 
degrees of freedom that can be set to zero by means of local 
unitaries in the first place, and by properly taking into account 
the symmetries of the problem. The optimization problem for 
the conditional entropy, and equivalently for the discord, is re- 
cast into a compact form that implies an elegant relationship 
among the eigenvalues of the ensemble obtained after the local 
measurement process on one qubit. The derived transcenden- 
tal constraints are amenable to direct numerical solution. 

We have then compared quantum discord with an alterna- 
tive but affine quantity, the geometric discord Dq, identifying 
the classes of states with extremal values of geometric discord 
at fixed quantum discord. For a fixed geometric discord, max- 
imal quantum discord is attained by different families of states 
depending on the degree of non-classical correlations, encom- 
passing pure as well as mixed, symmetric and nonsymmetric 
states. In general, the hierarchical bound D < y/lDc holds 
for all two-qubit states. 

We hope that our results can provide further insight into 
the fascinating but still not well-understood paradigm of 
non-classical correlations beyond entanglement in composite 
quantum systems. The methods presented here can be gen- 
eralized to higher-dimensional and continuous variable sys- 
tems. 
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